Decoherence of one-dimensional flying qubits due to their cross-talk and imperfections 
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We study decoherence of propagating spin- 1/2 excitations in generic (non-integrable and/or 
disordered) spin chains. We find the relevant decoherence times to be shorter in both the near- 
critical and diffusive regimes (if any), which fact might have important implications for the recently 
proposed spin chain-based implementations of quantum information processing. 



In spite of its success in describing thermodynamics 
and, more recently, spin/thermal transport, the theory 
of quantum one-dimensional {ID) spin systems has been 
struggUng to develop a systematic approach to such non- 
equilibrium problems as, e.g., the time evolution of a 
certain (in general, non-stationary) input state even in 
the simplest exactly solvable cases (see, e.g., Refs. [1]). 
However, it is exactly this kind of problems that have 
been brought about by the recent advents of spintronics 
and quantum information processing (QIP). 

For example, many theoretical discussions of an ab- 
stract multi-qubit (N S> 1) register invoke a generic spin- 
1/2 Hamiltonian [2] 
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while the usually proposed architectures have the topol- 
ogy of a ID, 2D or 3D array with the nearest-neighbor 
(NN) inter-qubit couplings (despite a desireable reduc- 
tion in the required number of qubits, the higher con- 
nectivity networks may not be that easy to assemble in 
practice) . 

In the QIP context, the issue of paramount importance 
is that of quantum coherence which, if no precautions are 
taken, can be easily destroyed by any coupling to an ex- 
ternal dissipative environment. Moreover, even though in 
the absence of such a coupling the density matrix p{t) of 
the entire A^-qubit system undergoes unitary evolution, 
any reduced density matrix (RDM) obtained by tracing 
out some of the degrees of freedom can still exhibit an 
ostensible decay caused solely by the qubits' "cross-talk" . 

In the QlP-related applications, one is primarily in- 
terested in the short-time behavior which strongly di- 
minishes any chances of the Poincare recurrences of the 
initial quantum state that might occur in the noiseless 
limit at longer times. Furthermore, an arbitrarily weak 
coupling to a noisy environment eliminates any possible 
quantum states' revivals and makes this "deterministic 
decoherence" a truly irreversible behavior. 

Contrary to the non-interacting case, however, the 
analysis of decoherence in a system of many coupled 
qubits described by Eq.(l) can no longer be carried out in 
terms of the conventional transverse/longitudinal relax- 
ation times (T1.2) of individual on-site spins. Instead, 
an adequate discussion must start out by identifying 



the proper basis of delocalized spin- 1/2 states with def- 
inite momenta (hereafter referred to as " spinous" ) that 
span the Hilbert space of the interacting multi-qubit sys- 
tem (1) and then proceed towards computing various n- 
spinon (n = 1,2,...) RDMs given by the partial traces 

T?'l,...,JV-n/5(t). 

Besides providing an appropriate framework for the 
discussion of decoherence, the spinon basis can also offer 
a natural implementation of the notion of "flying qubits" 
proposed as a vehicle for performing such important QIP 
tasks as teleportation or quantum key distribution [2]. 

To this end, in the present paper we evaluate the 
spinon decoherence rates in several situations that might 
be of interest for QIP and ascertain the prospects of using 
spinous in the QIP protocols relying on the availability 
of flying qubits. 

At first sight, the presence of unbound spinous in the 
spectrum of the Hamiltonian (1) may seem like a com- 
monplace, as suggested by the possibility to readily cast 
Eq.(l) in the form of the lattice Hamiltonian of spin- 
1/2 fermions with some short-ranged quartic interactions 
by applying a fermion representation of the spin op- 
erators {Si — ^ipl^Uaa'tpia') Complemented by a local 
single-occupancy constraint {tpl^ipi^ = 1). However, the 
answer to the question as to whether or not the "con- 
stituent" fermion operators ^ia correspond to any (pos- 
sibly, approximate) eigenstates of the Hamiltonian (1) 
would strongly depend on the actual form of the latter. 

The existence of deconfined spinous in higher dimen- 
sions (or a lack thereof) had been the subject of a long 
ongoing debate which was rekindled by the discovery of 
cuprates representing a new class of strongly correlated 
Mott insulators. Despite an extensive work on the topic, 
however, only a handful of solid arguments in favor of 
the occurrence of such excitations in a number of rather 
exotic models have been made so far (see, e.g., [3] and 
references therein). 

In contrast, many ID spin systems do conform to the 
picture of a "spinon Fermi liquid" (SFL) whose proper- 
ties appear to be quite different from those of the spin- 
ordered states described in terms of the conventional 
spin-1 magnons. Incidentally, the interest in spin-chains 
as a promising layout of the practical quantum register 
has been bolstered by the reports [4] of a massive buildup 
of quantum entanglement at those points in the parame- 
ter space of Eq.(l) where the system undergoes quantum 
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phase transitions. 

It should be noted, however, that, depending on the 
qubits' physical makeup, the corresponding Eq.(l) may 
lack any particular spin symmetry, thus hampering the 
applicability of such customary ID techniques as the 
Bethe ansatz, bosonization, and conformal field theory. 

As an alternate approach that has been previously em- 
ployed in various strongly correlated fermion systems, we 
adopt the semi-phenomenological spin- fermion model de- 
scribed by the action [5] 



k,q 



clEduj 



(2) 



written in terms of the continuum fermionic variables 
which are no longer subject to the local single-occupancy 
constraint, while their interactions are mediated by 

spin density fluctuations governed by the dynamical 
spin susceptibility x"^('^,a) = Y^ij I dte"^*~"'^'~'^ < 



The gapful fermion dispersion efc = Vv'^k'^ + ac- 
counts for the main effect of the processes of backward 
scattering, and the gap A <C J is assumed to vanish at 
certain values of the microscopic parameters in Eq.(l). 

Thus, the model (2) describes a gas of (for any A 0, 
massive) spinons that are still subject to a forward scat- 
tering, the singular part of which is controlled by the 
effective coupling constant g measuring a deviation from 
the parent SFL state. In turn, the latter is thought to 
be composed of quasiparticles that can be adiabatically 
related to the one-spinon states \k(j >= ipl.^\0 > even 
in the presence of non-singular residual interactions that 
were omitted in Eq.(2). 

Right at the quantum critical point (QCP) A = 0, and 
at all energies below a high-energy cutoff set by the av- 
erage exchange coupling J the T = susceptibility has a 
generic power-law form 



r]''''{q)J' 



2c(„ 



(3) 



In the absence of spin-rotational invariance, the expo- 
nents aab and the prefactors r}°'^{q) depend on the direc- 
tion in spin space. In the isotropic {SU{2) — ) or axially 
{U{1)-) symmetric case, however, all (or some) of the 
numerators vanish at zero transferred momentum q in 
compliance with the conservation of the corresponding 
component of the total spin. 

At temperatures above the gap (A < T < J) the 
susceptibility demonstrates a characteristic "quantum- 
critical" (QC) behavior 



q^e - - ilQcYI^^ 



(4) 



where w/^ - - 7qc cx T and xt{T) cx {J/Tf'^-^'- [6]. 

Notably, the temperature-dependent relaxation rate 
IQG = i-^TO^ lnx(w,0)|a;^o manifested by Eq.(4) can 
be present even in the host SFL state (i.e., for g = 0). 

At low temperatures (T < A) the off-critical suscepti- 
bility may also differ between the " quantum disordered" 
(QD) and "renormalized classical" (RC) regimes [6]. In 
either case, however, it can be cast in the form of a sum 
over the intermediate (1, 2, . . .-fermion) states, the coher- 
ent (one-particle) component Xco/j(c^,<z) assumes a form 
similar to Xqci'^^l) where Xo*'(^) is now given by the 
same expression as in Eq.(4) with the temperature substi- 
tuted by A, while the relaxation rate jqd.rc oz Te^^f^ 
[6,7] is suppressed as compared to its QC counterpart. 

In contrast to the pole-like structure of xtohi'^^l)^ ^ 
generic incoherent two-particle contribution 
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(5) 



features either a singularity {Jiah < 1) or a cusp (/?afe > 1) 
at the Lo = 2A threshold. 

In what follows, we consider a representative example 
of the spin-isotropic system supporting massive spinons 
that is provided by the frustrated Heisenberg model 
(J° = J,5j- ,±1 + J'<5,,,±2, Bf = 0). 

In this model (which can also be mapped onto a two- 
chain zigzag ladder), a sufficiently strong (J' > w 
0.24J) ncxt-nearest-neighbor (NNN) exchange gives rise 
to a spontaneous dimerization and opens up a gap A rx 
J exp{—const.J/ J' — J'^). 

For all J' < the system remains in the same criti- 
cal state as the standard NN Heisenberg model described 
by the SU{2)k=\ Wess=Zumino-Witten theory. Accord- 
ingly, the susceptibility given by the hydrodynamical ex- 
pression Xcoh('^'9) = 5'^^v'^q^ /{v'^q^ — w^) shows neither 
a gap, nor an anomalous exponent. 

The effect of the NNN and other frustrating ex- 
change interactions which render spinons massive can 
also be studied in the framework of the sine-Gordon (SG) 
model where the spinons can be identified with the SG 
(anti)solitons [8]. 

The incipient SFL instability with respect to dimer- 
ization occurs when the (anti)solitons become degener- 
ate with some of their bound states ("breathers"), thus 
resulting in the formation of massive triplets (prototypes 
of "flying qutrits"). 

In the gapped regime, the low-T susceptibility takes 
the form (5) where riab{q) oc q^ and = |(5afe, 

thereby indicating that the forward scattering remains 
non-singular and the (anti)solitons largely preserve their 
integrity. This expectation is supported by our calcula- 
tion of the spinon decoherence rates (see Eq.(12) below). 

Nonetheless, by analogy with the stability criteria of 
the conventional Fermi liquid, one might expect that a 
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singular (long-ranged and/or retarded) spin fluctuation- 
mediated interaction could make the spinous ill-defined 
and give rise to a critical "spinon non- Fermi liquid" state 
characterized by anomalous exponents. 

The subsequent analysis shows that such a behavior 
can indeed occur in spin- anisotropic systems, a popu- 
lar exactly solvable example of which is presented by 
the XF-model in transverse field (J? = — ^^(1 + 
S)Sj,i±t, Jl = - Ji(l - ^)5i,i±i, Bl = B). 

After the Jordan- Wigncr transformation followed by 
the Bogoliubov rotation, the XF-model with 5^0 gets 
transformed into a free gas of Majorana fermions with 
the velocity v = J6 and the gap A = B — J that vanishes 
at the Ising QCP [6]. 

These excitations possess no conserved physical spin 
and can be thought of as domain walls (Z2-kinks) be- 
tween two locally degenerate ground states. The fact 
that the kinks lack a simple representation in terms of 
the constituent fermions ipka implies that the correspond- 
ing one-spinon RDM p{x — y,t) =< ^l{x,t)'tpa' {y,t) > 
should exhibit strong decoherence (see Eq.(13) below). 

In support of this expectation, the critical susceptibil- 
ity of the model in question is now given by Eq.(3) where 
?7"^(0) 7^ and the anomalous exponents depend on 6 
(e.g., a^x = 3/4 and 7/8 for 5 = and 1, respectively). 

In the QD regime {B > J), the functions ^'^^'VV'^y 
feature a pole-like structure, while in the RC domain 
{B < J) they exhibit a kinematic square-root singularity 
{Pxx,yy,xy = 1/2) at the 2 A threshold. In contrast, 
displays a cusp (Pzz =3/2) and remains incoherent for 
all values of the transverse field [9] . 

The untying framework of the spin-fermion model (2) 
also allows one to investigate the effects of such qubits' 
imperfections as inhomogcncous local field or exchange 
couplings that vary from one site to another. In the pres- 
ence of a static disorder, the SU (2)-invariant low-energy 
{w,vq < 7(r) susceptibility acquires a diffusive pole 



(6) 



where the spin diffusion coefficient Dg = v^/2^i;r is de- 
termined by the disorder- induced transport rate jtr- It 
is worth mentioning that in a generic (non-integrable) 
spin chain Dg is finite at T > even in the absence of 
disorder, consistent with a non- vanishing spin (thermal) 
resistivity of such systems [10]. 

Similar to the critical spin fluctuations, the soft diffu- 
sive mode can prodive for a possible source of singular 
forward scattering between the spinons (see Eq.(14) be- 
low). 

In either case, the effect of the low-energy spin- 
fluctuation exchange on the spinon evolution operator 
U^{x,t) can be studied in the framework of the contin- 
uum equation 

[idt - e{idx) + iCT"S'"(a;, t)]U'^{x, t\S) = id{x)5{t) (7) 



where S{x, t) represents an instantaneous configuration 
of the effective spin environment that the propagating 
spinon is subjected to. 

The formal solution to Eq.(7) can be cast in the form 



JJ{x,t\S) 



dr 
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(T"S"'(w,g)(e'«^-l)] 



where the matrix ordering P is to be performed accord- 
ing to the time order in which the a-matrices appear in 
the series expansion of the exponential operator. 

Using Eq.(8) one can construct a one-spinon RDM 
given by the usual double-time Kcldysh path integral 
p{x - y,t\S) = JdudvU'^ix - u,t\S)p^-"^iu - v)U^{v - 
yT—t\S) for a given realization of the spin environment 
and an initial density matrix [x) . 

Upon statistically averaging f){x — y,t\S) over S{x,t), 
we arrive at the expression 



p{x,t) = J 



dyL{x-y,t)p'-''>{y)e 



(9) 



where we singled out the (super)operator of unitary evo- 
lution (hereafter z = x — y + iO) 

L{z,t) = 1^6*'=^ cos^ Cfei = 

k 



Z ztvt 

z^vt 



(10) 



At A = the kernel L{z, t) — 1 ^_|_ 5{z ± vt) describes 
a mere ballistic spreading of the entanglement present in 
the initial state (a "non-relativistic" counterpart of this 
behavior was studied in the case of the ferromagnetic 
Heisenberg model in Ref. [11]). 

Apart from the propagation, entanglement undergoes 
an apparent decay characterized by the decoherence fac- 
tor 



<^{z,t) 



2-K 



(w - ek, + ek,+q) 



/mx""(w, g)(coth ^ - tanh ^^^)(1 - cosqz) (11) 

where k^, « 1/z, and the spinon Fermi distribution fmic- 
tion eliminates contributions of the bosonic modes with 
energies u)>T which, otherwise, would have resulted in 
a residual decoherence at T = 0. 

Notably, Eqs.(9-ll) bear a certain resemblance to the 
results of the semiclassical calculation of Refs. [10] where 
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the one-particle RDMs of the exactly solvable trans- 
verse field Ising model and 0(3) non-linear cr-model were 
shown to factorize into a product of a quantum propa- 
gator and a relaxation factor, the latter becoming trivial 
at T = 0. 

Now, estimating Eq.(ll) on the semiclassical trajec- 
tory X — y = vt, wc can deduce a characteristic decoher- 
ence rate T from the equation F = -^^{vt,t)\t^i/r- 

In the spin-isotropic case, we obtain 

Tsym <x g^T, T > A 

which is comparable to or smaller than the SFL relax- 
ation rates jqc,QD,rc caused by a non-singular scatter- 
ing neglected in Eq. (2) , thus suggesting that the spinous 
remain well-defined elementary excitations. 

By contrast, in the totally anisotropic case the deco- 
herence rate is strongly enhanced as compared to Eq.(12) 

Tasym CX g2 j2aj.l-2a > ^ 

^S^y^e-^/-, T<A (13) 

The decoherence rates given by Eq.(13) become greater 
than the SFL ones at max[T,A] < cx .g^/"J which 
signals a rapid decay of any spinon-likc wave packet with 
a definite spin S = 1/2, thus indicating the difficulty 
of creating robust flying qubits with the spin- anisotropic 
Hamiltonians. 

At temperatures below T* the analysis of decoherence 
requires one to proceed beyond the lowest order of the 
cluster expansion employed in the derivation of Eqs.(9- 
11). 

Besides, Eqs.(12,13) demonstrate that the decoherence 
exhibited by the one-spinon RDM can be suppressed 
by the spectral gap and, conversely, it tends to become 
stronger as the system is tuned towards a QCP. This con- 
clusion complements the earlier analysis of the decoher- 
ing effect of an external dissipative environment which, 
unlike in the present case, remains unchanged [12]. 

The decoherence rate also appears to be enhanced in 
the disordered case where the susceptibility takes the 
form (6) 

F.., oc g'/'T'/Hl' (14) 

The diffusive contribution (14) dominates over its bal- 
listic counterpart (12) at all temperatures up to T*, cx 
Itr/g'^ where both terms are of order jtr, while below 
T*** oc g'^'jtr the rate r^if becomes of order T and a sys- 
tematic account of higher orders of the cluster expansion 
is again required. 

Besides being manifested by the single-spinon RDM, 
the rate (14) controls various amplitudes that are sen- 
sitive to the spinon phases such as, e.g., the "spin- 
Cooperon" amplitude associated with the pairs of time- 
reversed trajectories that return to the starting point. 



Thus, T(iif can also be extracted from the weak- 
localization correction to the low-temperature spin (ther- 
mal) conductivity 




^y/2^-Hvt,t)^^ oc (i).rd,/)V2 (15) 



that should be contrasted with the lowest order (Drude) 
result ai^^ oc D^Xu- 

It was argued in Refs. [13] that in the QIP context the 
qubits' localization by disorder may prove advantageous, 
as far as the goal of creating well-localized spin- 1/2 states 
is concerned. However, our results show that increasing 
the amount of disorder would come at a price of greater 
decoherence rates. 

To summarize, in the present work we studied the ap- 
parent spin decoherence that occurs at T > in even a 
closed non-integrable ID spin system. This short-time 
"deterministic decoherence" can be revealed by, e.g., the 
RDMs of a finite number of spinons propagating in a 
dissipative environment created by the rest of the inter- 
acting qubits. 

We focused on the near-critical and diffusive regimes 
where the emergence of soft collective modes (critical and 
diffusive, respectively) is more likely to cause a singular 
forward spinon scattering and an ensuing rapid decay of 
the spinon RDMs. 

In the case of a spin-isotropic Eq.(l) and no disorder, 
we found that spinon remains as good of a coherent ele- 
mentary excitation as it is in the parent SFL model. By 
contrast, in the spin-anisotropic case and/or in the pres- 
ence of spin diffusion the decoherence rates are strongly 
enhanced, consistent with the fact that the system does 
not support coherent spin- 1/2 states with definite mo- 
menta. We also observed that the decoherence rates can 
be reduced by tuning the system away from an incipient 
criticality and/or localization. 

Lastly, our findings suggest that the enhanced deco- 
herence (as quantified by the spinon RDMs) in the near- 
critical and diffusive regimes might significantly reduce 
any possible benefits of both the emergent massive en- 
tanglement [4] and the disorder-induced localization of 
the flying qubits' [13] that have been recently discussed 
as potentially attractive features of the spin chain-based 
QIP designs. 
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